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ABSTRACT 

We present a catalog of over 130,000 quasars candidates with NIR photometric properties, with an 
areal coverage of approximately 1,200 deg 2 . This is achieved by matching the Sloan Digital Sky Survey 
(SDSS) in the optical ugriz bands, to the UKIRT Infrared Digital Sky Survey (UKIDSS) Large Area 
Survey (LAS) in the near-infrared YJHK bands. We match the w 1 million SDSS DR6 Photometric 
Quasar catalog to Data Release 3 of the UKIDSS LAS (ULAS), and produce a catalog with 130,827 
objects with detections in one or more NIR bands, of which 74,351 objects have optical and if-band 
detections and 42,133 objects have the full 9-band photometry. The majority (~ 85%) of the SDSS 
objects were not matched simply because there were not covered by the ULAS. The positional standard 
deviation of the SDSS Quasar to ULAS matches is <5r.a. = 0.1370" and S Dccl . = 0.1314". We find an 
absolute systematic astrometric offset between the SDSS Quasar catalog and the UKIDSS LAS, of 
|R. A. ff sc t | = 0.025", and |Decl. offset | = 0.040"; we suggest the nature of this offset to be due to the 
matching of catalog, rather than image, level data. Our matched catalog has a surface density of « 53 
dcg~ 2 for K < 18.27 objects; tests using our matched catalog, along with data from the UKIDSS DXS, 
implies that our limiting magnitude is i ~ 20.6. Color-redshift diagrams, for the optical and NIR, 
show the close agreement between our matched catalog and recent quasar color models at redshift 
z < 2.0, while at higher redshifts, the models generally appear to be bluer than the mean observed 
quasar colors. The gJK and giK color-spaces are used to examine methods of differentiating between 
stars and (mid-redshift) quasars, key to currently ongoing quasar surveys. Finally, we report on the 
NIR photometric properties of high, z > 4.6, and very high, z > 5.7, redshift previously discovered 
quasars. 

Subject headings: catalogs - quasars: general 



1. INTRODUCTION 

With the completion of the Sloan Digital Sky Survey 
(SDSS; I York et al.lkOOOD . the era of high quality, homo- 
geneous CCD imaging over large fractions of the sky has 
arrived. The final da ta release from the SDSS, (DR7; 
lAbazaiian et al.ll2009D contains over 11,000 deg 2 of imag- 
ing data, with 357 million unique objects being identified 
in a ~ 60 Terabyte database. 

The identification of quasars was a major focus of 
the SDS S project; utilizing the broad 5-filter pho- 
tometry (IRichards et al.l 120021) led to efficient selection 
(jVanden Berk et all 120051: IRichards et al.l I2006D of low, 
z < 2.2, and high, z > 3.5, reds hift quasars, identified 
via their spectroscopic signatures (jSchneider et al.ll2010L 
and references therein). However, quasars can also be 
efficien tly identified by the ir SDSS imaging properties 
alone ([Richards et alj 12001). due to their point-source 
appearance, but non-stellar locat ion in color-color space . 
With the most recent catalog of IRichards et al.1 ([2009a, 
hereafter R09), based on imaging from the SDSS, the 
number of photometrically identified quasars now stands 
at over 1 million. Both the SDSS spectroscopic and 
photometric quasar catalogs have been used to inves- 
tigate global qua sar properties such as the luminosity 
function (QLF; iFan et al.l [20011 IRichards et all [2001 
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| Croqm et all l2009aD and clustering (IMyers et al.l 120061 
20071 iShen et al.ll2007l . 120091: iRoss et al.ll2009f ). 

The emergence of large surveys has not been confined 
to the optical regime. In the near-infrared (NIR: A ~ 
1 — 5/im) quasar catalogs fe.g. IBarkhouse fc Hal|[2001 : 
I Cutri et all [2001 iFrancis et all 12004 IQfek et all 12007: 
IKouzuma fc Yamaokall2010al lbf) have been available since 
the completion of the 2 Micron All Sky Survey (2MASS; 
ISkrutskie et a l. 2006). Table Q] summarizes how previous 
quasar surveys, both in the optical and NIR, compare by 
area, number of objects and magnitude range. 

Quasar observations in the NIR are particularly im- 
portant for individual objects that are seen only in 
the reddest, or indeed potentially none, of the opti- 
cal filters (e.g. IFan et al.l 120061: iVenemans et al.l 120071: 
IWillott et al.ll2009D . Observations of the general quasar 
population in the if-band are key since this links the 
rest-frame ultraviolet (UV)/optical to the mid-infrared 
(MIR; A « 5 — 30/xm); the former being where there is 
the peak in the radiative output for Type I , non-obscured 
quasa r s, (e.g. iShakura fc Sunyaevlll973t iSanders et al.l 
[19891: iKishimoto et al.l 120081: IRichards et al.1 12009b). 
and the latter where reprocessed light heats intrinsic 
dust to ~ 30 - 300 K (e.g fPier fc Krolifl li"99l 
lEfstathiou fc Rowan-Robinsonl [19951: lLacv et al.l I2004D . 
Observations in the observed JC-band also measures the 
rest-frame i and g-bands at redshifts z ~ 1.9 and z ~ 3.7, 
respectively. 

However, it is only the bright, g < 16 magnitude 
quasars that are detected by the relatively shallow limits 
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Table 1 

Descriptions of previous catalogs of quasars. a O/I, Optical or Infrared; P/S, Photometric or spectroscopic survey. ^SDSS Faint Quasar 
Survey. c The FIRST-2MASS Red Quasar Survey. d The SDSS DR7 Quasar catalog provides NIR photometry from 2MASS for 53,584 
objects, 29,551 or which are detected in the if-band. e The 2MASS Second Incremental Data Release. 



Survey 


Area (deg^ ) 


Nn 


Magnitude Range 


O/I & P/S a 


Reference 


COMBO-17 


0.8 


192 


R < 24 


O/P 


Wolf et al. ( 20031 


SFQS 6 


4 


414 


g < 22.5 


o/s 


Jiang et al. (20061 


SDSS-ULAS DR1 


189 


2 873 


K < 18.2 


I/S 


Chiu et al. (20071 


2SLAQ QSO 


190 


8 764 


18.00 < g < 21.85 


o/s 


Croom et al. (2009b) 


2QZ 


700 


23 338 


18.25 < bj < 20.85 


o/s 


Croom et al. (20041 


SDSS-ULAS DR3 


1 200 


74 351 


K < 18.4 


I/P/S 


this work 


FIRST-2MASS RQS C 


2 716 


57 


K < 14.3, (R - K ) > 4, ( J - K) > 1.7 


I/S 


Glikman et al. (20071 


SDSS DR6pQ 


8 342 


1 015 082 


i < 21.3 


O/P 


Richards et al. (2009a1 


SDSS DR7Q d 


9 380 


29 551 


K < 17.0 


I/S 


Schneider et al. 2010 


SDSS DR7Q 


9 380 


105 807 


i < 20.2 


o/s 


Schneider et al. 2010 


2MASS 2IDR e 


~ 30 000 


2 277 


K < 15 


I/S 


Barkhousc & Hall (2001) 



of the 2MASS survey, and the majority of known quasars 
are fainter than this in the NIR band s. The UKIRT 
Infra red Deep Sky Survey (UKIDSS; ILawrence et al.l 
I2007I) . a seven-year sky survey which began in 2005 
May, has five different survey components. The UKIDSS 
"Large Area Survey" (ULAS) aims to reach ~ 4 magni- 
tudes deeper than 2MASS over an area of up to 4,000 
deg 2 , directly overlapping the optical imaging footprint 
of the SDSS. 

The primary goal of this paper is to cre ate a catalog of 
over 130,000 quasars with optical, ugriz (|Fukugita et al.l 
1996), and NIR, YJHK photometry, with an areal cov- 
erage of 1 200 deg 2 . By matching the catalogs of R09 in 
the optical regime to that of the ULAS in the NIR we 
produce, by a factor of at least two, the largest catalog 
of quasars detected in the if-band. The major motiva- 
tion for the catalog will be its utilization in an future 
study where we measure the iC-band quasar luminosity 
function (Peth, Ross et al. in prep.). By concentrating 
on the w 200 deg 2 area from the SDSS known as "Stripe 
82", our i*T-band quasar luminosity function will show 
the evolution from redshift of zero to two. 

The analysis by iTrammell et al.l (|2007j) is an exam- 
ple of the synergy produced by the matching and pro- 
duction of a multi-wavelength catalog. These authors 
match ^6000 SDSS quasars to UV data provided by 
the G alaxy Evolution Explorer (GALEX; IMartin et al.l 
2005) satellite and find that over 80% of the optically de- 
tected quasars have near-UV detections. The quasars are 
well separated from stars in UV-optical color-color space. 
The large sample size allows for the construction of SEDs 
in bins of redshift and luminosity, which shows the me- 
dian SED becoming bluer at UV wavelengths for qu asars 
with lower continuum luminosity. IBall et ail (|2007[ ) also 
perform catalog matching using SDSS quasar data and 
GALEX UV photometry, with the goal of understanding 
quasar photometric redshift prope r ties. 

The st udies by I Warren et all (|2000H . iCroom et ail 

(120011). iSharo et all (120021) and m ore recently, 

iMaddox et al.l (|2008l ) and iSmail et all (|2008l ). are 
another motivation why a sample of quasars with 
NIR photometric properties is desired. These authors 
show that using a "KX selection", where the quasar 
SED shows an excess in the iT-band compared to a 
stellar SED, can successfully identify quasar candidate 



objects that would be normally excluded from the SDSS 
(optical) quasar selection algorithm - even for dust 
reddened quasars. This is an important result since 
selecting complete quasar samples via the KX method 
opens up the possibility of investigating the "Quasar 
Epoch" over the redshift range of 2.2 < z < 3.5, where 
current, usually optically selected, quasar samples are 
particularly p o orly r epresen ted. A sim il ar pro ject to 
IMaddox et all (|2008l ) was iNakos et all (|2009l ). who 
also select quasar candidates using the KX-technique, 
where these authors identify quasars on the basis of 
their optical (R and z ) to NIR (K«) ph otometry and 
point-like morphology. lJurek et al.l (|2008l ) also test the 
KX method and find that it is more effective than the 
traditional "UV Excess" (UVX) selection method at 
finding red, (bj — K) > 3.5, quasars. 

There ar e two comparable stud ies and samples to ou r 
own work. iChiu et all (|2007l ) and lSouchav et all (|2009f ). 
the latter recently producing the "Large Quasar Astro- 
metric Catalog" (LQAC). Our work differs from these 
studies, in two key ways; (i) we have over an order 
of magnitude mo re objects in our sample compared to 
IChiu et all (|2007l ); (it) our catalog uses UKIDS S data, 
as opposed to 2MASS data (Souchay et al. 20Q9|), where 
the former is much better matched to the SDSS imaging 
depth. 

This paper is organized as follows. In Section 2 we 
present our sample, giving a brief overview of the SDSS 
and UKIDSS. Section 3 lays out our catalog. In Section 4 
we present N(z), color-redshift and color-color relations 
from our matched catalog. In Section 5 we present anal- 
ysis of high-redshift quasars and calculations of the i and 
if-band number counts. The Appendix gives further de- 
tails on magnitude conversions and cross-checks of our 
study. 

We assume the currently preferred flat, "Lambda Cold 
Dark Matter" (ACD M) cosmology where ftb =0.042, Q, m 
= 0.2 37, fl A = 0.763 ([Sanchez et all 120061 : ISpergel et all 
12007( 1 and quote distances in units of hr 1 Mpc to aid 
in ease of comparisons with previous results in the liter- 
ature. Where a value of Hubble's Constant is assumed, 
e.g. for absolute magnitudes, this will be quoted explic- 
itly. All optical magnit udes are based and quoted on the 
AB zero-point system (|Oke fc Gunnll 1983ft . while all NIR 
magnitudes are based in the Vega system, with conver- 
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sions from AB to Vega given in the Appendix. 

2. DATA AND METHODS 

In this section, we provide overviews of the surveys and 
the catalogs we utilize, and how matching was performed. 

2.1. The SDSS DR6 Photometric Quasar Catalog 

Details regarding the SPSS can be found in the s eries 
of SDSS Data Release papers ( Abazaiian et al.ll2009l and 
references therein). Full details for the SDSS DR6 Photo- 
metric Quasar Catalo g (DR6pQ) are given in R09, with 
IRichards et al.1 (120041) providing the "proof-of-concept" 
study and lWeinstein et all (|2004h providing an empirical 
algorithm for obtaining photometric redshifts. Here we 
present the details specific to our work. 

The photometric imaging data for the DR6pQ 
is based upon the SPSS Data Release 6 
(|Adelman-McCarthv et all 12008ft . Points sources 
with PSF i-band magnitudes between 14.5 and 21.3 
are extracted from the SPSS Catalog Archive Server 
(CAS). We continue the conven tion of R09, utiliz- 
ing iibercalibrated magnitudes (|Padmanabhan et al.l 
2008) which are available in the SPSS database. The 
iibercalibrated magnitudes represent the most robust 
photometric measurements as they are calibrated across 
SPSS "stripes" to a single uniform photometric system 
for the entire SPSS area. All magnitu des have been 
correc ted for Galactic extinction using the lSchlegel et all 
(1998) dust maps. 

There are 1 015 082 objects in total in the R09 catalog 
across 8 342 deg 2 . The PR6 primary imaging data cover 
an area of 8417 deg 2 , although as noted in R09, due to 
cuts the total effective area covered by this catalog is 
reduced by ^75 deg 2 . 

2.2. The UKIRT Infrared Deep Sky Survey 

lLawrence et al.l (|2007[ ) gives the general overview for 
UKIPSS. In brief, the UKIPSS is a collection of five 
survey s of different cover age and depth and using WF- 
CAM (|Casali et al.ll2007h on UKIRT. WFCAM has an 
instantaneous field of view of 0.21 deg 2 , and the various 
surveys employ up to five filters, ZYJHK, covering the 
wavelength range 0.83-2.37^m. The photome t ric sy stem 
and calibration are de scribed in lHewett et al.l (|2006t ) and 
iHodgkin et al.l (|2009( ). respectively. The pipeline pro- 
cessing is described in Irwin et al. (2010, in prep.) and 
the W FCAM Science Archive (WSA) by lHamblv et all 
(2008). The processed right ascension and declination 
data are accurate to 0.1 arcsec. We have used data from 
the worldwide 4th data release, PR3, which is described 
in detail by Warren et al. (2010, in prep.). For this pa- 
per, we concentrate on the Large Area Survey and the 
Peep Extragalactic Survey. 

2.2.1. UKIDSS LAS 

The ULAS aims to map a large fraction of the North- 
ern Sky, ~ 4000 deg 2 , which, when combined with the 
SPSS, produces an atlas covering almost an octave in 
wavelength. The target depths of the survey are Y = 
20.3, J = 19.5, H = 18.6, K = 18.2 (Vega), and the 
ULAS does not image in the WFCAM Z-band. The 
ULAS data for our matched catalog came courtesy of 
private communications with Mike Read from the WSA. 



Unlike the SPSS, the ULAS multiband pho tometry is 
not taken in one observation (e.g. Sec. 5.2 of iPve et ahl 
I2006t lLawrence et al.1 120071 Sec. 4.2). Therefore the 
four bands can, and do, have different coverage maps, 
with the H and K bands obtained together, and Y and 
J done separately. The ULAS PR3 coverage is 903 deg 2 , 
1,161 deg 2 , 1,091 deg 2 and 1,111 deg 2 in Y, J, H and K 
respectively. Note that the ultimate 7- year ULAS goal is 
to cover 4028 deg 2 , in each filter, and have two passes of 
the entire ULAS area with the J filter. The PR3 cover- 
age of each of the four bands, and for all four bands, is 
shown in Figure [TJ The joint coverage of all four bands 
is 801 deg 2 . 

2.2.2. UKIDSS DXS 
The UKIPSS Peep Extragalactic Survey (PXS) plans 
to cover 35 deg 2 of sky to a 5er point-source sensitiv- 
ity of J=22.3 and ^=20.8 in four specifically selected 
multiwavclcngth survey areas. The locations and areas 
of these four fields, XMM-LSS, Lockman Hole, ELAIS 
N l and SA22 (a. k .a VI MOS-4), are given in Table 5 
of lLawrence et ahl (|2007h . Three of these, the Lockman 
Hole, ELAIS Nl and SA22 fields, are covered by the R09 
PR6pQ catalog. We shall therefore use the PXS data to 
test the faint end and limiting magnitude of the optical 
and NIR matched quasar catalogs. 

2.3. SDSS Stripe 82 

Stripe 82 is a 300 deg 2 area of repeat photometry on a 
2.5 degree wide stripe, centered on the celestial equator in 
the Southern Galac tic Cap and running fro m 300° to 60° 
in R.A. (Sec. 3.2, lAbazaiian et al.ll2009D . Co-addition 
of the best of this data means that optical imaging on 
Stripe 82 can reach roughly 2 magnitudes fainter than 
the main survey and as we s hall see in section 15.21 this 
deeper imaging was used by Uiang et "all (|2009l ) to dis- 
cover new, very high, z ~ 6, redshift quasars. We do not 
require this deeper imaging for the preparation of our 
SPSS-ULAS matched catalog. However, due to the fact 
that there is almost complete coverage in the ULAS, as 
well as many other multi-wavelength surveys, along with 
the regular rectangular geometry (that will simplify the 
number counts and luminosity function calculations in 
future work), we utilize heavily the data from this area. 

2.4. Matching Procedure 

Our matching procedure is conceptually straight- 
forward, though is made relatively cumbersome by the 
shear number of objects to process. Before the million 
strong catalog (PR6pQ) was uploaded, tests were per- 
form using the "CrossIP" form. The positions of the 
1,015,082 objects from the PR6pQ catalog were uploaded 
to the WFCAM Science Archive. To count as a match, 
the returned NIR object must be the closest object in the 
ULAS to the SPSS coordinates, selected from within a 
circle of radius 1" (though we note a more stringent 0.5"' 
matching radius returns essentially identical results) . As 
long as at least 1 NIR band contained a non-default, i.e. 
not -9.9999xl0 8 value, the object was considered a good 
match. Table [5] shows the total number of objects in the 
R09 catalog and the number of those objects that are 
matched to ULAS and PXS separately, both for the full 
PR6pQ coverage and specifically Stripe 82. 
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Figure 1. Coverage over the entire UKIDSS area in each of the 4 NIR bands and finally in all 4 NIR bands simultaneously. Uses no 
projection, with axes in degrees for both Right Ascension and Declination. Displays YJHK coverage in descending order, followed by 
coverage in all four bands. 
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Figure 2. The differences in Right Ascension and Declination 
recorded by SDSS against that recorded by the ULAS for all 17,835 
matched objects in Stripe 82. A circle of radius 1" is shown. 



Table 2 

Overall Catalog Numbers. 



Catalogs 


Full DR6 


Stripe 82 


Ri09 


1,015,082 


36,625 


DR6pQ-ULAS Matches (in >1 NIR band) 


130,827 


17,835 



Table 3 

Spurious match counts. Results of secondary matching tests 
where we offset the right ascension and declination of Stripe 82 
DR6pQ catalog objects and then upload them to WSA to observe 
the probability of spurious matches. 



Offset 


R.A. offset "matches" 


Decl. offset "matches" 


1" 


7253 


10965 


2" 


56 


40 


5" 


75 


72 


10" 


84 


93 


60" 


76 


78 



We also perform matching tests where we offset the 
right ascension and declination of Stripe 82 DR6pQ cata- 
log objects and then upload them to WSA to observe the 
probability of spurious matches. Table |3] shows the re- 
sults of our matching tests. Once the offsets are increased 
to more than 1" the amount of spurious "matches" drops 
off precipitously. 

Figure [2] displays the separation in arcseconds between 
our matched objects in the SDSS and UK1DSS. The his- 
tograms, with a bin size of O.Of " running along the sides 
of the plot, represent the distribution of the Right As- 
cension and Declination separations. We find the stan- 
dard deviation of these histograms to be Sr. a. = 0.1370" 
and ^Doci. = 0.1314". These values compare very well 
to those dete rmined in the U KIDSS Early Data Releas e 
paper (EDR; iDve et al.|[2006h and bv lChiu et all (|2007[) . 



Also similar to lChiu et al.l (|2007f) . we find an absolute off- 
set between the SDSS and ULAS of |R.A. offsct | = 0.025", 
and |Decl. or i S ot | = 0.040". We suggest the origin of this 
offset to be due to the matching of catalogs, rather than 
image level data, and is not due to the astrometric cali- 
bration of either of the surveys, nor the size of CCD pixels 
(the SDSS camera pixel size is 24 /im/0".396 on the sky; 
WFCAM pixels are 18/im/0".4 on the sky). However, we 
also suggest this issue warrants further investigation. 

3. THE MATCHED CATALOG 

Our matched catalog between the SDSS DR6pQ and 
the ULAS DR3 (hereafter simply referred to as "the 
matched catalog") comprises a total of 130,827 objects 
matched between the R09 catalog and the ULAS. The 
first ten objects are given in Table |4l with the column 
titles, meanings and format given in the following text 
and in Table [5j The full version is found in the electronic 
version in machine readable form. 

Columns 1 through 4 The first four columns des- 
ignate the internal catalog number; the formal name of 
the quasar as reported from the SDSS, with the form 
at SDSS Jhhmmss.siddmmss.ss; and the ObjID num- 
bers specific to the SDSS and UKIDSS, specifically from 
the PhotoObj All and lasSource database tables respec- 
tively. 

Columns 5 and 6 give the J2000 right ascension and 
declination from SDSS in degrees. 

Columns 7 through 9 shows the relative positional 
measurement accuracy between SDSS and UKIDSS. 
AR.A. and ADecl. values are calculated by subtracting 
positional measurement as recorded by SDSS by the po- 
sitional measurement as recorded by UKIDSS. Distance 
values are not calculated explicitly, rather they are val- 
ues reported directly from the WFCAM database. The 
"distance", i.e. difference, is simply given by: 

d = \J (R.A. sdss - R.A.tjkidss) 2 + (DccI.sdss - Decl.uKiDSs) 2 

Distances are the only angular quantities expressed in 
arcseconds rather than degrees. 

Columns 10 through 19 give the PSF magnitude 
values in each of the 5 SDSS optical bands, and their 
associated error. Errors are not re-calculated, but are 
those previously reported in R09. 

Columns 20 through 27 provide NIR magnitudes, 
jAAperMagS, yAperMagS, hAperMagS and kAperMagS, 
and the associated error as given in the ULAS database. 
The AperMag3 magnitudes are the aperture corrected 
magnitudes measured by UKIDSS, which records all the 
flux in a 2.0" diameter and that for typical seeing con- 
ditions, pro vide the most ac curate estimate of the total 
magnitude (|Dve et al.1 12006). The jAAperMagS magni- 
tude is reported, since the ULAS ultimately aims to have 
J-band imaging in two epochs. For quasars lacking cover- 
age in one or more bands, the default reported magnitude 
and error are -9.9999 x 10 s . Table[6]shows the breakdown 
of matches based on each NIR band and combination of 
bins. Objects numbers for Stripe 82 only are given in 
parentheses. 

Columns 28, 29, 30 and 31 give the nominal 5-cr 
depth, uid, at the given coordinates for a R09 quasar cat- 
alog member that falls within the ULAS DR3 footprint 
(in any one of the four NIR bands) . Following equation 
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Table 4 

The first ten objects in our NIR-matchcd quasar catalog. The catalog will be published in its entirety in the electronic edition of the 

Astrophysical Journal. The first 10 rows, of 130,827 total, are shown here. 



No. 


Name 


SDSS ID 


ULAS ID 


R.A.QnQQ 


Decl.cncic 


A R a. 




Diff 


1 


000001.38-010852.2 


588015507658768592 


433825841152 


0.0057816 


-1.1478427 


-0.180 


0.001 


0.180 


2 


000001.93-001427.5 


588015508732510344 


433822400512 


0.0080576 


-0.2409745 


-0.130 


-0.137 


0.189 


3 


000006.42+005206.3 


588015510343123177 


433816469504 


0.0267757 


0.8684207 


0.031 


-0.316 


0.318 


4 


000006.53+003055.2 


588015509806252150 


433817550848 


0.0272281 


0.5153489 


0.118 


0.164 


0.202 


5 


000007.58+002943.3 


588015509806252166 


433817550848 


0.0316036 


0.4953744 


-0.039 


0.127 


0.133 


6 


000008.13+001634.6 


588015507658768713 


433819123712 


0.0338984 


0.2763040 


-0.077 


0.108 


0.132 


7 


000009.31-010703.1 


588015507658768768 


433825644544 


0.0387928 


-1.1175526 


-0.050 


0.105 


0.117 


8 


000011.96+000225.3 


588015509269381139 


433820860416 


0.0498422 


0.0403718 


-0.064 


-0.078 


0.101 


9 


000012.25-003220.5 


587731185667080338 


433823776768 


0.0510761 


-0.5390492 


0.033 


-0.316 


0.317 


10 


000012.27-010405.5 


588015507658768732 


433825841152 


0.0511651 


-1.0682041 


-0.062 


0.193 


0.203 



Table 5 

Format of our Matched catalog. With field names and descriptions. 
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a 
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A Decl 


F12 7 


L-'lllCICilCC J3CIWCCI1 L-'CCl. 1IUII1 OiJOO dllU UlVlJJOO ^diCbeCb^ 


9 


Diff 


F12.7 


Difference between TJKIDSS Object and SDSS Object designated as a Match (arcsecs) 


10 


u 


F12.7 


u SDSS PSF magnitude in AB 


11 


u err 


F12.7 


u Error SDSS PSF magnitude in AB 


12 


9 


F12.7 


g SDSS PSF magnitude in AB 


13 


g err 


F12.7 


g Error SDSS PSF magnitude in AB 


14 


r 


F12.7 


r SDSS PSF magnitude in AB 


15 


r err 


F12.7 


r Error SDSS PSF magnitude in AB 


16 


i 


F12.7 


i SDSS PSF magnitude in AB 


17 


i err 


F12.7 


i Error SDSS PSF magnitude in AB 


18 


z 


F12.7 


z SDSS PSF magnitude in AB 


19 


z err 


F12.7 


z Error SDSS PSF magnitude in AB 


20 


Y 


F12.7 


Y ULAS PSF magnitude in Vega 


21 


Y err 


F12.7 


Y Error ULAS PSF magnitude in Vega 


22 


J 


F12.7 


J ULAS PSF magnitude in Vega 


23 


J err 


F12.7 


J Error ULAS PSF magnitude in Vega 


24 


H 


F12.7 


H ULAS PSF magnitude in Vega 


25 


H err 


F12.7 


H Error ULAS PSF magnitude in Vega 


26 


K 


F12.7 


K ULAS PSF magnitude in Vega 


27 


K err 


F12.7 


K Error ULAS PSF magnitude in Vega 


28 


depth.y 


F12.7 


nominal 5— a depth in the Y-band 


29 


depth_J 


F12.7 


nominal 5— <r depth in the J-band 


30 


depth_ff 


F12.7 


nominal 5— a depth in the /f-band 


31 


depth.K 


F12.7 


nominal 5— cr depth in the ii"-band 


32 


zPhot 


F12.7 


Photometric Redshift 


33 


zPhotLow 


F12.7 


Lower limit of the photometric redshift range 


34 


zPhotHi 


F12.7 


Upper limit of the photometric redshift range 


35 


zSpec 


F12.7 


Spectroscopic Redshift (when known, —1 otherwise) 


36 


zphotprob 


F12.7 


Probability Photometric Redshift is accurate 


37 


E(B-V) 


F12.7 


Galactic Extinction 



Table 6 

Counts for each of the NIR Bands in the catalog, along with colors and count in all 4 NIR Bands in bold. The italic numbers in 

parentheses are the totals for Stripe 82 only. 



Band 


Y J 


H 


K 


J + H J + K 


J + H + K 


Y+ 
J+ 
H+ 
K+ 


80,544(13,681) 59,085(9,384) 
89,962 (11,141) 


54,965 (8,690) 
51,701 (8,236) 
72,347 (11,191) 


53,656 (8,442) 
50,264 (7,993) 
61,015 (8,892) 
74,351 (11,546) 


47,787 (7,258) 46,338 (6,938) 


42,133 (6,007; 
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(5) from lDve et al . (2006), md is the 5cr detection limit 
for a point source, denned as five times the standard de- 
viation of the counts in an aperture, corrected for flux 
outside the aperture: 

m d = to - 2.51og 10 (5cr(1.2A'p) 1/2 /£ OX p) - TO ap 

where too is the photometric zero-point (the PhotZPCat 
attribute in the WSA); <7 s k y is the measure of the stan- 
dard deviation of counts in the sky (skyNoise); N p is 
the number of pixels in the aperture, the factor of 1.2 ac- 
counts for the covariance between pixels and t e xp is the 
Exposure Time, (expTime). m ap is the aperture correc- 
tion term, AperCor3, which gives the correction needed 
for a 2.0" aperture diameter (i.e. the respective quantity 
for the quoted aperMag3 values). These values can all be 
found in the Multif rame and Multif rameDetector ta- 
ble of the WS A. Section 1 0, and in specific section 10.2.2, 
example 3, in lDve et al.l (|2006|) has further details here. 

Columns 32 through 34 represent the photomet- 
ric redshif10 fr om R09, with columns 29 and 30 giving 
the lower and upper limits of the photometric rcdshift 
range respectively. These values should be taken with the 
zphotprob quantity, described below. Note that these 
upper and lower redshift ranges a nearly, but ultimately 
not, symmetric about the reported photometric redshift, 
and therefore should not be treated as the formal error 
on the photometric redshift. 

Column 35 gives the spectroscopic redshifts are from 
the R09 catalog, and are ultima tely based on matche s 
to the SDSS DR5 quasar catalog dSchneider et aT1l2007l ). 
the 2QZ quasar catalog ()Croom et al.l 120041 ) . the 2dF- 
SDSS LRG and QS O Survey (2SLAQ) quasar catalog 
(ICroom et alJl2009bD. and the SPSS DR6 spectroscopic 
database (jAdelman-McCarthv et al.l I2008D . By design 
the R09 catalog does not have spectroscopic information 
for all of its objects. We find that out of the 130,827 
objects matched from R09 and ULAS, 20,740 have spec- 
troscopic redshifts as reported in R09, 6,623 of which are 
on Stripe 82. The value —1 is given when there is no 
corresponding spectroscopic redshift. 

Columns 36 represent the probability of the reported 
photometric redshift being in the given redshift range, 
with R09 (their section 4.6 and Figure 15), and sec- 
tion 4 in our investigations, giving further details on the 
zphotprob quantity. 

Column 37 represents the estimated Galactic red den- 
ing at the given position from ISchlegel et ail (|1998| ). 

The smaller UKIDSS DXS DR3 matched catalog, with 
1,070 matched objects, has exactly the same format as 
given here for the SDSS-ULAS matched catalog given 
here, and again will be available with the electronic ver- 
sion. The only notable difference being that since the 
DXS does not cover the F-band, all those magnitudes 
are reported as 0.00. 

4. GLOBAL CATALOG PROPERTIES 

3 On occasion, a photometric redshift of -1 is returned in the 
catalog; this is simply an error saying that the photo- z code failed 
for some, usually tractable reason (G.T. Richards, priv. comm.). 
This issue only affects 3,574 objects out of the total 1,015,082 from 
R09 quasar catalog and only 242 quasars on Stripe 82. We do not 
attempt to deal with this issue, leaving the resolution for the next 
version of the photometric catalog. 
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Figure 3. Distributions of each of the four NIR bands from Stripe 
82 ULAS matches in Vega magnitudes, Y {solid black), J, (dotted 
red), H (dashed blue) and K (dot-dashed purple), for the matched 
objects, normalized at the peak to unity. 

4.1. Coverage 

Only around 10-14% of the R09 objects match to the 
ULAS DR3, with the primary reason for this low percent- 
age being the difference in coverage between the SDSS 
DR6 and the ULAS DR3 (Figured]) gives the ULAS DR3 
coverage). To investigate how many R09 objects actu- 
ally lie with- in the ULAS DR3 footprint, we take the R09 
catalog and the "multiframe" information connected to 
UKIRT via the WSA. A multiframe is essentially the 
footprint of the four detectors of an individual WFCAM 
exposure (see e.g. section 10, iDve et aLll2006| ). In total 
there were 172,186 R09, optically detected objects that 
fall within the ULAS DR3 footprint, in the union of the 4 
bands. For individual bands, there are 109,959, 141,304, 
132,788 and 135,257 matches in the Y, J,H and if bands 
respectively, with 97,541 R09 objects were in the inter- 
section of the Y, J, H and K coverage footprint. Thus, 
and using Tabled we find that 73%, 64%, 54% and 55% 
of the SDSS R09 objects are detected in the DR3 foot- 
print, in the Y, J, H and K bands respectively. 43% of 
the objects have detections in the all four NIR filters. 

The surface density is «122 deg~ 2 for the R09 DR6pQ. 
This compares to the surface density for objects with 
K < 18.27 (see discussions below regarding this limiting 
magnitude) of «80 deg ~ 2 . These va l ues ar e generally in 
line with those given in iSmail et al.l (|2008l ) who derive a 
surface density of QSOs with K < 20 of between 85-150 
deg -2 . 

4.2. Magnitude Distributions 

Figure [3] shows the normalized magnitude distributions 
for the matched objects in each of the four NIR bands, 
showing at which magnitude each particular band is lim- 
ited by. As expected, each band peaks at a different 
value with shorter wavelength bands peaking at fainter 
magnitudes. 

For Figure [H the solid (red) histogram represents the 
objects which are matched in any of the 4 NIR bands, 
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Figure 4. Histogram showing population by i-band PSF magni- 
tude. Each line was normalized to a peak value of 1.0 and a bin 
size of 0.1. The purple dotted line shows the histogram based on 
objects from R09 in Stripe 82, the solid red line displays the his- 
togram of only ULAS matched objects and the black dashed line 
represents the histogram of DXS matches over the entire sky. 

while the dotted (purple) histogram represents all the 
optically detected objects from the DR6pQ on Stripe 82 
only. The dashed (black) line are the matches to the 
deeper DXS, the limiting; depth of which is K = 20. 8 
(|Lawrence et al.ll2007t lStottll2007t iSwinbank et~aLll2007[) . 
Due to the smaller amount of available data for UKIDSS 
DXS, the matches characterized by the dashed line is less 
smooth than the other two histograms. The sharp cut-off 
at i=21.3 displays the i-band limit in the original SDSS 
DR6 photometric redshift catalog. 

The i-band magnitude distribution of the matched cat- 
alog strongly suggests that the ULAS is not complete at 
the optically faint end of DR6pQ. The suggested com- 
pleteness limit for the matched catalog is thus i w 20. G. 
We also caution that our matched catalog is a heteroge- 
nous sample of objects, and not a complete statistical 
sample as is. This will have particular ramifications when 
the calculation of the NIR Quasar Luminosity Function 
(QLF) comes to be calculated in our companion paper 
(Peth, Ross et al., in prep.). The significant shift be- 
tween the ULAS matched histograms and the other two 
(optical-Stripe 82 and DXS matched) distributions most 
likely implies brighter optical objects are preferentially 
matched to ULAS NIR detections. It is reassuring to 
see that the DXS and optical-only histograms seem to 
match very well in shape and limiting magnitude, and 
the matched DXS catalog is provided as a test-bed for 
future investigations into the fainter, i > 20.6 end of the 
.fT-band quasar population. 

Figure [5] gives the color-magnitude distributions for 
both the i (black points) and K (red points) bands. 
Again the sharp cut-off at i = 21.3 is seen. The 
.ftf-band limit is seen to be down to «18.4, poten- 
tially a little deeper, consistent with Figure [3] This 
is slightly deeper than th e preliminary depths reported 
in lLawrence et al.l (|2007ft but in good agreement with 



Figure 5. Color-magnitude histograms for both the iab an d the 
^*Vega bands, for i — K color. Note the lack of points with faint 
(i > 20) blue, (i — K) < 2.5, colors and the effect of the limiting 
A'-band magnitude for redder, (i — K) > 2.7, objects. 

the up dated (UKIDSS DR2) values from iWarren et al.l 
(2007) and indeed is also consistent with the mea- 
surements from the ULAS team, who give a 5cr point 
source limiting magnitude of 18.27 (R.G. McMahon, 
priv. comm.). We suggest these faintest if -band ob- 
jects are actually 5<r detections in a non-K NIR band, 
but then also have some flux in the aperture used to cal- 
culate the kAperSMag, just possibly not to as strong a 
limit as the quoted 5er for point sources. Also apparent 
are the clear selection affects associated with the matched 
catalog, with few faint i-band objects at (i — K) < 2 and 
very few objects fainter than K = 18 and redder than 
(i — K) > 3. Again, the location of the matched cata- 
log objects in e.g. K vs. (i — K) color-magnitude space, 
will become critical for statistical calculations using these 
data. 

As a check, we visually examine (via the SDSS CAS 
Image List Tool and the UKIDSS WSA Getlmage form) 
the very reddest objects in our catalog. These very red 
objects are defined as having if-band detections and col- 
ors of (r— K) > 5.0 or (i—K) > 4.5. There criteria return 
17 and 10 objects respectively. Upon visual inspection of 
these 17 (10) sources, 6 (4) turned out to be either asso- 
ciated with, or contaminated by a foreground extended 
source. However, 13 out of the 17 (r — K) selected ob- 
jects had suspiciously high, z = 4.5 — 5.6, redshifts, or 
zphotprob < 0.8 values, indicating the deceptive nature 
of these objects. Of the 10 (i — K) selected sample, the 4 
sources associated with the foreground galaxies, also had 
zphotprob < 0.8 values, while the remaining 6 objects ap- 
peared as genuine point sources. These particular very 
red (i—K) objects deserve future follow-up investigations 
and we shall discuss in more detail in the next section 
the utility of the zphotprob flag. 

4.3. Redshift Distributions 

Figure [6] shows the redshift distribution, N(z), for 
quasars, in redshift bins of Sz = 0.1. Three scenarios 
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Figure 6. Distributions of Photometric and Spectroscopic Redshifts, for the full DR6 coverage (upper left); Stripe 82 only (upper right); 
iC-band matches over the DR6 (lower left). The black line represents Photometric redshift distributions and the orange line represents 
spectroscopic redshift distributions. The subsample of quasars from Stripe 82 (the upper-right plot) posses nearly identical redshift 
distributions to those over the entire sky. However, objects with i<"-band coverage have photometric redshifts that are preferentially in the 
low redshift bins. The difference in numbers given here and in Table [6] are due to the sample number of objects which have a photometric 
redshift of zPhot= — 1. 



are presented; the top left panel is for the full, optical- 
only, SDSS DR6 photometric redshift catalog; the lower 
left panel is from the matched catalog with i^-band de- 
tections across the full coverage, and top right is from 
the matched catalog with X-band detections from Stripe 
82. In each panel, the redshift distribution of the photo- 
metric redshifts is given by the solid (black) line, while 
those objects with spectroscopic redshifts are shown by 
the dashed (orange) histogram. 

For the optical-only distributions the redshift distribu- 
tions are potentially well matched, apart from the excess 
of photometric z ~ 2.3 objects which distort the normal- 
ization. 

The redshift histograms for those objects with if-band 
detections from the matched catalogs, for both the full 
coverage and Stripe 82, are almost identical. Interest- 
ingly, the spike in the number of objects at z ~ 2.3 



seems to be dramatically reduced (although does not dis- 
appear completely) for the AT-band matches. As we shall 
see in § 4.6, when we examine the stellar contamination 
of our sample using the giK color-color space, we find 
that (i) stellar contamination is potentially higher in ob- 
jects with photometric redshifts z ~ 2.3 and (ii) due 
to the respective shape of typical quasar and (e . g. M - 
type) stellar SEDs (e.g. Figure 1, iMaddox et aLll2008l) , 
quasars will be preferentially selected over stars, if mea- 
sured in the if-band. Thus, we suggest the z ~ 2.3 spike 
is because these mid-redshift objects have a lower "zphot- 
prob" value than the general photometric quasar sample. 
As explained later in § 4.6 we see that the stellar locus 
separates nicely for z < 3.2 but not for higher redshifts. 

Once the -ft"-band matches are made, the photometric 
and spectroscopic histograms now seem to be in reason- 
able agreement, though there are possibly decrements 
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Figure 7. Left panel: The spectroscopic redshifts plotted against photometric redshift for all matches with spectroscopic redshifts. Black 
points have no zphotprob cuts while red points have zphotprob> 0.8. The 1:1 relation is plotted with a purple line. Right panel: Histograms 
of difference between zspec and zphot for (1) NIR matches over the entire sky (upper right); (2) NIR matches on Stripe 82 (middle right); (3) 
K- matched objects on Stripe 82 (lower right). A small full width half maximum says that the majority of photometric and spectroscopic 
redshifts are in agreement. 



of photometric objects at z ~ 0.7 and more severely 
at z ~ 1.8. We note that this is not a new feature, 
having been seen in R09 (their Figure 14), and as such 
are motivated to continue investigations for the photo- 
metric and spectroscopic redshift distributions from our 
matched catalog. 

Figure [7] shows the relationship between the photo- 
metric redshifts and the spectroscopic redshifts for the 
objects in our matched catalog. From top to bottom, 
the left panel of figure [7] is; the entire matched cat- 
alogue with spectroscopic coverage; all NIR matches 
on Stripe 82 and all if -band matches on Stripe 82. 
A purple line shows the 1:1 relation between spectro- 
scopic and photometric redshifts. There is considerable 
structure in the photometric-spectroscopic quasar red- 
shift relation, which ag a in has been seen previously, e.g. 
IMountrichas fc Shanks! IpOOH ). R09. Dc generacies be- 
tween z ~ 0.8 and z ~ 2.4, as well as the width of 



the distributions around z ~ 1 — 1.5 can be understood 
by the appearance, and potential confusion of, various 
stron g emission lines in the various (SP SS) filters (see 
e.g. rF^r][i999t iRichards et all [20021 l200i for detailed 
descriptions). There are a couple of notable differences 
between the left hand panels of Figure [7] that should be 
briefly highlighted. First, note that a number of objects 
with high, z > 3, photometric redshifts and no zphotprob 
cut, disappear when there is a match to the NIR data. 
Upon checking, we suggest this is caused by these objects 
generally being lower luminosity AGN than "standard" 
quasars, and thus thought to be at higher- z than was 
spectroscopically found, and in general just too faint for 
the ULAS to detect. Also, though, this is good motiva- 
tion for the utilization of the zphotprob flag. A second 
feature is the continued structure in the NIR matched 
plots. This, coupled with the general lack of structure 
at z > 3, suggests that as far as photometric redshift 
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Figure 8. Color versus redshift plots from Stripe 82 matches of u — g, g — r, r — i and i — z. Mean magnitude (red line) in a binsize of 
0,1 mag, 1(7 dispersion (striped yellow line) both above and below the mean values, and Model Quasar colors (orange line) as determined 
m IHewett et al.ll2uM) are each plotted in terms of redshift. 



estimation for quasars (and not e.g. differentiation be- 
tween high- z quasars and cool stars), then the directly 
matched NIR photometry, i.e. "Is there simply NIR pho- 
tometry?" , actually does not provide a great amount of 
extra information. Of course this does not take into ac- 
count the actual detected NIR fluxes (or upper limits) 
that would help refine photo- z estimates, or the use the 
NIR to find very high, z > 6, redshift objects. 

The right panel of figure [7] shows the histogram and 
gaussian fits with calculated full width half maxima for 
the difference in photometric and spectroscopic redshifts. 
We define that quasars with photometric and spectro- 
scopic redshifts that have a difference greater than So- 
to be "catastrophic failures". All quasars with a differ- 
ence in photometric and spectroscopic redshifts that fall 
within the 3er are deemed "non-catastrophic failures" . 
The full numbers on the agreement between spectro- 
scopic and photometric redshifts can be seen in table 

m 

Red points in the left panel of figure [7] show quasars af- 



ter a redshift probability, zphotprob cut of 0.8, the value 
of 0.80 being inspired by the checks done by R09 (their 
Figure 15). This cut reduces the amount of catastrophic 
failure redshifts by a considerable margin, although the 
FWHM of the (zSpec-zphot) difference only changes very 
marginally. Although the zphotprob> 0.80 is a "blunt" 
cut, we suggest this as a relatively sensible division, but 
do note ther e are r ecent, more sophisticated studies, e.g. 
IMvers et all ()2009D . that show how one can take into 
account the full photometric redshift probability distri- 
bution function when performing statistical calculations. 

4.4. Optical and NIR Colors 

We now turn to the optical and NIR colors of our 
matched catalog. 

Figure [5] shows the optical colors, (u — g), (g — r), (r — i) 
and (i — z), as a function of redshift, for our matched 
catalog on Stripe 82. The mean colors of the data are 
given by the solid (red) line, and were determined in 
a given bin of width 5z — 0.1, with the associated la 
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Table 7 

The number of objects with photometric and spectroscopic that agree to within a 3<r difference. The total number of objects in each bin 

is displayed in italics. 





SDSS DR6 


% 


Stripe 82 


% 


K-matched 


% 


zprob > 0.0 


15527 (20605) 


75 


4933 (6589) 


75 


4218 (5644) 


75 


zprob > 0.8 


8981 (10188) 


88 


2403 (2620) 


92 


2114 (2292) 


92 



standard deviations given by the striped (yellow) line. 
The orange triangular-pointed line represent the "Aver - 
age" QSO model co lors presented in lHewett et al.l (|2006f ). 
IHewett et all (j2006f ) also present "blue" and "red" model 
quasar spectra, where the power-law slope, a, at wave- 
lengths A <12,000A is a = 0.0,-0.6 and -0.3 for 
blue, red and average model quasars respectively. We 
cut the plots at redshift z — 3, since we feel we have 
poorer statistics at higher redshifts and that there are 
not enough points to render a representative mean value. 
We will examine the colors of individual high-z quasars 

in Section [5J 

Pre vi ous studies fe.g.lRichards et al.ll200lUJester et all 
2001 iRichards et all [200l ICroom et al.l I2009bl: 
Hewett fc Wild! 1201 OH have provided detailed stud- 



ies of the features of quasar colors in the SDSS color 
space, and their form as a function of redshift. As 
such, we will not repeated those analyses here. Instead, 
we will make comp arisons to the models presented in 
Hewett et all (120061) of which details can be found in 
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Starting with (u — g) , there are potentially two redshift 
ranges, z = — 0.6 and z > 2.3 where the mean colors of 
the matched quasars are generally redder than the model 
colors. While the differences are relatively small, ~ 0.2 
mag at the low redshift end, ~ 0.1 mag at the high red- 
shifts, the model tracks are consistently bluer than the 
observed quasars, and at z < 0.6, barely consistent with 
the la spread of the quasar (u — g) color. Note also 
however, that care has to be taken in the determination 
of the (u — g) color at z > 2.2 since this is where the 
lyman-a forest first enters the optical bands, and the 
determination of accurate it-band photometry becomes 
increasingly problematic. 

For (g — r), the model tracks reproduce the mean 
quasar colors well at all redshifts, z — — 3. However, 
there might be a tend for the observed quasars to be gen- 
erally redder at low (z < 0.5) and high (z > 2.3) redshift, 
though this can be seen as marginal. 

For (r — i), it appears that the model tracks do in 
fact reproduce the mean quasar colors very well, for z < 
2, and especially for z = 0.7 — 2.0. The disagreement 
between model and observed mean at z « 0.5 seems to 
be caused by the fact that there are a noticeable number 
of (red and blue) outliers, with the observed colors of 
the general population better tracked. As can be seen, 
the model colors agree well with the raw data at z « 
0.5, (r — i) ~ —0.1. The model tracks begin to struggle 
to reproduce the observed colors at z > 2, but remain 
consistent with the la dispersion until z > 2.5. 

For (i — z), just as for (r — i), the model tracks replicate 
the mean colors very closely, all the way to z ~ 2.7 and 
thus for the optical colors, the models perform best for 
these reddest bands. At z > 2.7 for (i — z) and indeed 
all four optical colors, the models struggle to reproduce 
the observed quasar colors at these higher redshifts, con- 



sistently being too blue. 

We tentatively suggest t his is the same affect a s has 
recently been reported in [P rochas ka et all (|2009| ) and 
iWorseck fc Prochaskal (|2010l ). namely, that the original 
SDSS color-selection that was used to select quasars, and 
from which the R09 and hence our own matched cata- 
log is based, may be systematically biased as these red- 
shifts. It has long been known that the redshifts around 
z ~ 2.7 are very troublesome for selecting quasars, since 
the broadband colors are so simi lar to thos e of e.g. F5 
V stars (see e.g. Figure 1 of iFanl I1999T) . leading to 
very poor survey effi c iency and heavy in completeness 
(|Richards et all 120061: ICroom et all I2009bl ). However, 
these recently analyses model in great detail the spectra 
of these high redshift objects, paying particular atten- 
tion to the flux transmission (and decrement) blueward 
of Lyman-a, in order to generate accurate mock UV and 
optica l photometry. Specifically, IWorseck fe Prochaskal 
(2010) suggest that the SDSS color-selection selection 
systematically misses quasars with blue, u — g < 2 col- 
ors at 3 < z < 3.5 due to the preference of selecting 
z > 2.7 quasars with intervening H I lyman limit sys- 
tems (which generally turn quasars red) . This would po- 
tentially begin to explain the offset between the model 
tracks and observed quasar colors in Figure [5] We are 
keen to note however, that there are more recent v e rsions 
of the quasar color models given in lHewett et all (2006) 
and that further analysis is required here. 

Noting the shape, response and wavelength coverage 
(and indeed gaps between filters, see e.g. Figure 1 
IHewett et al.l l2006h between the Y, J, H and K filters, 
we can now begin to explain the NIR color trends with 
redshift seen in Figure [U The NIR color-redshift rela- 
tions still exhibit structure, but arguably less than in 
the optical. This is mainly due to the fact that it is 
now the lines towards the redder end of the optical, and 
in particular Ha, t hat lead to the tr ends seen in Fig- 
ure [5] _Wejiote that | Assef et all (|2010f ) and in particular 
Glik man et al.l (|2006l their table 6) are useful references 
for the following. 

The Ha line is very much the dominant emission line 
over the rest-wavelength range, A = 6500 — 21000A. The 
blended He I and Paschen-7 lines at around 10,900A, 
the Paschen-/? line at 12,820A, and Paschen-a emission 
at 18,756A, are also present but contribute little. For 
example, even if a Paschen line had an equivalent width 
of 100A, with the NIR bands being so broad (e.g. ~ 
4000A for the if-band; see Table [5]) , the contribution 
here would be ~a few % at redshift 2 = 0- significantly 
smaller than the measurement errors in the colors. These 
emission lines are superimposed on the continuum power- 
law spectrum, F(v) oc v a where a = —0.78, measured 
over A = 5700 - 10850A ([Glikman et al.| [2fM but see 
also lKishimoto et all (120081 )1 

Thus, starting with the top left panel, (z — Y), we note 
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Figure 9. Color versus redshift plots from Stripe 82 matches, z — Y (top left); Y — J (top right); J — H (bottom left) and H — K (bottom 
right). The given lines are the same as in Figure [51 Note the dip and rise at z ~ 0.4, ~ 0.7, ~ 0.9 and ~ 1.3, for z ~ Y, Y — J, J — H and 
H — K respectively, due to the location of the strong Ha emission line in the NIR filters. 



the dip at redshift z ~ 0.3 is due to Ha entering the z- 
band causing the general quasar color to be blue, while 
at redshift, z ~ 0.45, Ha is leaving z, and entering the 
Y-band causing the colors to become redder, reaching 
an observed mean (model) peak of (z — Y) « 0.5 (1.0) 
at redshift z w 0.6. At z ~ 0.7, Ha leaves the Y-band. 
The model tracks trace the observed colors well up to 
z ~ 2, but then grow progressively redder and redder at 
higher redshifts (though do remain within the measured 
la dispersion). This behavior, the model tracks being 
redder than the observed quasar colors at high, z > 2, 
redshift is only see in (z — Y). 

Moving to, (Y — J), we again see a dip, at z ~ 0.5, 
towards bluer colors, as Ha enters the Y-band, and then 
a rise towards redder colors, peaking at redshift z ~ 0.9 
as Ha moves out of Y, and into the J-band (noting the 
gap between bands at 1.4/im). The models' continued 
rise and fall behavior at z ~ 1.1 — 2.1 can be explained 
by H/3, then O III and then Mg II marching through 
the WFCAM bands. The observed mean colors however, 



show a much smaller amplitude of color-change over this 
redshift range, though the individual objects seem to re- 
spond more to these emission lines. 

For (J — H), the passing of Ha first through the J- 
band and then through the .ff-band gives the trough at 
z ~ 0.8 — 1.1. Here the model tracks, even out to our 
maximum plotted redshift of z = 3, replicate the ob- 
served quasar colors very well. The success of the model 
tracks in (J — H) could potentially be used for quasar 
selection in future redshift surveys, inclu ding the im- 
plementation of "KX" completeness tests (|Croom et al.l 
[20011: ISmail etalll200l Uurek et alJl2051 . 

Finally, for {H — K), the Ha trough is across redshifts 
z = 1.2 — 1.8. Note the broadening of the trough in 
the four NIR colors as Ha passes to higher redshifts. 
Again the model tracks replicate the observed quasar 
colors well, though there is a large dispersion at all red- 
shifts. The one place the models tracks are in disagree- 
ment is at low, z < 0.5, redshifts. This is potentially due 
to flux from an underlying red host galaxy being over- 
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represented, but more likely, the lack of data points to 
calculate the observed mean well. 
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Figure 10. The location of matched quasars, and stellar sources 
in the g — J vs. J — K plane. Stars are solid (black) points; 
< z < 1 quasars are open (blue) triangles; 1 < z < 2.2 quasars 
are solid (purple) circles; 2.2 < z < 3.2 quasars are open (orange) 
squares and (r ed) crosses are z > 3.2 objects. The ("average") 
model tracks of Hcwctt ct al. (2006) are given by the (green) line, 
with z = around (J — K) = 2.5 increasing in Az = 0.1 steps 
(green diamonds) to z = 4.8 at (g — J) = 8. 



4.5. Stellar Contamination 



Foll owing iMaddox et "all (|2008t ) and iSmail et al.1 
(2008), we now investigate the location in (g — J) vs. 
(J — K) ("gJK") color-space of our matched catalog. 

Here our motivation is not so much directed at creating 
complete "KX" selected samples, but is more towards de- 
termining if there is a certain part of color-space where 
quasars, and in particular mid, z — 2.2 — 3.5, rcdshift 
quasars lie apart from stars. This is of key importance 
to two new spectroscopic surveys, the SD SS-III: Baryon 
Oscil lation Spectroscopic Survey (BOSS; Schleg el et al.l 
120071) and the AAT-UKIDSS-SDSS (AUS) QSO Survey 
(P.I. S.M. Croom). Both the BOSS and AUS survey aim 
to gather spectroscopic information for z > 2.2 quasars, 
for investigations into cosmological parameters using the 
Lyman-a forest (BOSS) and global quasar population 
studies (AUS and BOSS). As such, any technique to re- 
duce stellar contamination is of upmost interest to these 
survey teams. The KX technique is based upon the fact 
that at redshifts z s»= 2 — 3, NIR photometry samples 
the Rayleigh-Jeans tails of A and F stars (the major con- 
taminants in an optical-based quasar selection), where 
flux is decreasing rapidly with the increasing wavelength, 
whereas quasar SEDs are remaining relatively flat. 

In Figure [TU] we plot the (g — J) vs. (J — K) col- 
ors for the objects on Stripe 82 from our match catalog, 
with spectroscopic redshifts. We also select 5000 point 



sources lying on Stripe 82, from the SDSS CAS, with 
spectroscopic redshifts less than 0.02, and gather ULAS 
NIR photometry for these objects where available. These 
can be considered stars because of their low redshifts and 
morphology, and are solid (black) points in gJK hgure. 
< z < 1 quasars are open (blue) triangles; 1 < z < 2.2 
quasars are solid (purple) circles; 2.2 < z < 3.2 quasars 
are open (orange) squares and (red) crosses are z > 3.2 
objects. Aga i n we show the "Average" model tracks of 
iHewett et atl(|2006l ). given by the (green) line, with z = 
around (J—K) = 2.5 increasing in Az = 0.1 steps (green 
diamonds) to z = 4.8 at (g — J) = 8. 
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Figure 11. g — i vs. i — K diagram to show the colors of Stripe 82 
objects from the matched catalog. Stars are solid (orange) points. 
Quasars are displayed in bins of varying redshift ranges: z < 1.0, 
black stars; 1.0 < z < 2.2 blue open triangles; 2.2 < z < 3.2 orange 
open squares, and greater than 3.2, red crosses. 

Our results are q ualitatively very similar to that of 
IMaddox et alj (|2008l ). that is, (a) the stellar locus traces 
out a clear band in gJK color-space, related to the dif- 
ferent stellar spectral type and (6) that in general, stars 
appear to have gJK colors that are in a distinct re- 
gion of color space than quasars, especially those with 
2.2 < z < 3.2 (orange) open squares in Figure [TU] We 
find this very encouraging and suggest that selection in 
the gJK color-plane could be of good utility to the BOSS 
and AUS surveys. We also note the almost complete 
lack of quasars in our gJK plot that have colors similar 
to those of (model) galaxies around {J — K) rs 1.7 and 
(g — J) m 4.5. However, a strong caveat that has to be 
imposed here is to recall that the original R09 quasar se- 
lection was predominantly geared towards z < 2 "UVX" 
selected, non-extend ed objects. T his is to be compared 
with the study by ISmail et all (|2008t l who find that 
lower, z r~j 0.45 compact galaxies are the main and ma- 
jor contami n ent o f their "KX" selection. Interestingly, 
ISmail et al.l (|2008f) report that the spectral mix of the 
contaminating population includes sources whose spectra 
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are best fit by both breadline and narrow-line AGN tem- 
plates, as well as narrow emission-line, absorption-line 
and spiral galaxy templates. Indeed, early analysis from 
the BOSS (Ross et al., 2010, in prep.) suggests that these 
morphologically compact, narrow-emission line galaxies 
at lower, z » 0.2 — 0.5 redshift, continue to be selected 

as z > 2.2 qu asars. 

Inspired bv lMehta et al.l ()2010[ ). we also plot the loca- 
tion of Stripe 82 objects from the matched catalog, this 
time in g — i vs. i — K (giK) color-space, see Figure ITTl 
Again, quasars are displayed in bins of varying redshift 
ranges: z < 1.0, black stars; 1.0 < z < 2.2 blue open tri- 
angles; 2.2 < z < 3.2 orange open squares, and greater 
than 3.2, red crosses. On initial inspection, we again see 
the stellar locus inhabiting a distinct sequence in color- 
space, which is somewhat separate to that of the quasars. 
However, we leave further investigations into the use of 
the (g — i) vs. (i — K) color-space for star-quasar-galaxy 
identification, and survey completeness to future study. 

5. HIGH REDSHIFT QUASARS 

5.1. High Redshift objects in Matched Catalog 

Figure [T^] plots high-redshift quasars, those with z > 
4.6 (essentially r-band drop-outs), in iYJ color-color 
space. Only objects with both a "good" photometric 
and spectroscopic redshifts are plotted. For a photomet- 
ric redshift to be considered "good" it must be within 
the error limits of the spectroscopic redshift. Out of 49 
high-redshift (z > 4.6) quasars possessing both a photo- 
metric and spectroscopic redshift, 40 have "good" pho- 
tometric redshifts. Out of this group of 40 quasars, 23 
objects possess both Y and J magnitudes, and we plot 
these in Figure [12] These objects span the redshift range 
of 4.606 < z spcc < 5.289. No obvious trends are visi- 
ble, arguably due to our small sample size. However, we 
can say that these high-redshift quasars have roughly the 
same color and foll ow the given high-z model tracks of 
IHewett et all (|2006| ) well. 

5.2. Very High Redshift objects 

For our final investigations, we actually leave the 
matched quasar catalog, and concentrate solely on the 
ULAS DR3 to investigate the NIR photometric prop- 
erties of the very high, z > 5.7, redshift objects that 
are given in th e series of Fan et al., and Jiang et al . 
papers papers, (|Fan et al.l 12001 IJiang et al.l 120081 . 120091 
respectively and references therein). We also note the 
19 very high redshift quasars discovered b y the Canada- 
France High-z Quasar Survey (CFHQS; IWillott et al] 
I21M I2TM l2Tjl0fK However, the CFHQS is not covered 
by the ULAS DR3, and only has the one (J) out of the 4 
ULAS bands, so the NIR properties of these objects are 
not reported here. 

Out the 30 objects - 19 from the Fan et al. studies, 11 
from Jiang et al. wit h one object in commo n and the ad- 
ditional object from IMortlock et al.l ()2009D . 13 (12) are 
within the ULAS DR3 (if -band) footprint, and 6 were 
detected in one or more of the ULAS bands. The photo- 
metric properties of these 6 detections, along with their 
spectroscopic redshifts, are given in Table [5] 5 of the 6 
objects have if-band detections, with magnitudes rang- 
ing from K = 16.19 — 18.16. The reason for only 5/12 
having if-band detections we suggest is primarily due to 



the limiting magnitude for the ULAS. As a quick check, 
we note that the Jiang et al. objects were discovered 
using the deeper SDSS optical imaging, this means, po- 
tentially 6 out of the 7 if -band non-detections are fainter 
than the ULAS if-band limit. The m ean J-band magni - 
tude from the 6 objects presented in IJiang et al.l (2009) 
is J « 21.0 ± 0.2, whereas for the objects in Table H] it 
is J rs 18.9 ± 0.3. Thus if the colors of these z ~ 6 
objects are (J — K) ~ 1.1 — 1 . 8, eve n the brightest and 
reddest object in IJiang et al.l ([2009) would be right on 
the ULAS if-band limit, with the average object being 
closer to K = 19.3. 

Our results are generally consistent with 
iVenemans et al.l (|2007l ) who report the discov- 
ery of the quasar at z = 5.85 (designated ULAS 
J020332.38+001229.2) afte r analysis of 106 deg 2 of sky 
from UKIDSS DR1, and IMortlock et all |2009) who 
report the di scovery of ULAS J1 319+0950, a quasar 
at z = 6.13. iGlikman etaD (|2008t ) surveyed 27.3 deg 2 
of the ULAS EDR, but found no z > 6 objects from 
spectroscopy of a 34 candidate list. 

Deeper and more complete surveys will be needed 
to get a better estimate on the spatial density of 
high-redshift quasars. To this affect, the Visible 
and Infrared S urvey Telescope for Astronomy (VI STA; 
lEmerson et"all 120041 : lEmerson fc Sutherland! 120101 ) has 
recently been commissioned, and is currently ramping u p 
its suite of six imaging surveys ([Arnaboldi et al.l 12010) . 
Of direct interest to the high redshift Universe will be the 
"VISTA Hemisphere Survey", (VHS; PI: R.G. McMa- 
hon), which aims to image the entire ~20 OOOdeg 2 South- 
ern Sky, down to Jab = 21.2 and K s A b = 20.0 (J « 20.3 
and K * 18.1). On a l onger timescale, t he Synoptic 
All-Sky Infrared (SASIR; iBloom et alll2009[ ) is planned. 
The hope and goals of these surveys will be to find 100's 
(VISTA) or potentially 1000's (SASIR) of z > 6 quasars. 
Whether these numbers are borne out in the new obser- 
vations, will, of course, be of great future interest. 

6. CONCLUSIONS 

We have matched the optical, photometrically selected 
quasar catalog of R09 to detections from the ULAS DR3. 
The positions of the 1,015,082 objects from the DR6pQ 
catalog were uploaded to the WFCAM Science Archive. 
To count as a match, the returned NIR object must have 
been the closest object in the ULAS to the SDSS coor- 
dinates, selected from within a pairing radius of 1". As 
long as at least 1 NIR band contained a non-default, i.e. 
not -9.9999 x 10 s value, the object was considered a good 
match. 

Our final catalog has 130,827 objects with detections in 
the optical and one or more NIR bands; of which 74,351 
objects have if-band detections and 42,133 objects have 
the full 9-band photometry. Using this matched catalog, 
we present the following main conclusions: 

• The positional standard deviation of the SDSS 
Quasar to ULAS matches is 5r.a. = 0.1370" and 
<!>Dcci. = 0.1314". We find an absolute systematic 
astrometric offset between the SDSS Quasar cat- 
alog and the ULAS, of |R.A. offsct | = 0.025", and 
Decl. ff sc t I = 0.040"; we suggest the nature of this 
offset to be due to the matching of catalogs, rather 
than image level data. 
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Table 8 

Photometric properties of 6 very high z > 5.7, redshift quasars as detected in the ULAS DR3, along with their spectroscopic redshifts. 
Positions are in degrees and J2000 coordinates. References: "Venemans et al. (2007). 6 Jiang et al. (2008). c Fan ct al. (2006). d Mortlock 

et al. (2009). 



R.A. 



Dccl. 



redshift 



Y 



J 



H 



K 



"• b 30.8849 0.2081 

6 58.45719 1.0680 

c 129.1827 0.9148 

d 199.7970 9.8476 

c 212.7970 12.2937 

c 245.8826 31.2001 



5.850±0.003 
6.049±0.004 

5.82±0.02 
6.127±0.004 

5.93±0.02 



26.21±0.57 25.80±0.55 22.75±0.23 



24.77±0.99 24.79±0.71 23.22±0.28 



23.72±0.22 20.87±0.10 19.85±0.12 
24.03±0.30 20.54±0.08 20.12±0.16 
21.23±0.09 18.86±0.05 18.26±0.03 
22.55±0.09 19.99±0.03 19.22±0.06 
22.86±0.30 19.65±0.08 19.47±0.07 



..22±0.02 23.96±0.91 25.44±0.62 24.88± 0.58 24.50±0.62 20.07±0.10 



19.06±0.10 
19.46±0.16 
17.70±0.029 
18.69±0.05 
19.22±0.093 
19.17±0.111 



17.76±0.07 17.32±0.09 
18.55±0.17 18.16±0.22 
17.03±0.03 16.19±0.03 
17.55±0.08 
18.29±0.09 17.46±0.08 
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Figure 12. High, 4.6 < z < 5.35 redshift quasars in i — YvsY — J 
color space. The "average" quasar model tracks of IHewett et al.l 
(2006) are again given, starting at z = 4.6 (labelled) and then 
increasing in Az = 0.1 steps, shown by the open (red) diamonds. 

• Our matched catalog has a surface density of ^108 
deg -2 , for objects detected in any of the four NIR 
bands, and w 53 deg -2 for objects K < 18.27. This 
compares to the « 122 deg~ 2 found from the R09 
DR6 pQ catalog and 85- 150 deg" 2 down to K < 20 
from lSmail etall (l200l) . 

• Tests using our matched catalog, along with data 
from the UKIDSS DXS, implies that our limiting 
magnitude is i « 20.6 and that the reddest, (r — 
K) > 5.0 objects turn out to be either associated 
with, or contaminated by, a foreground extended 
source. 

• We plot redshift histograms for the Stripe 82 sub- 
sample, i-T-matches sample and total matched sam- 
ple. The photometric N(z ~ 2.3) spike seen for the 
"total sky" sample appears to be dramatically re- 
duced (although does not disappear completely) for 
the JT-band matches. The photometric and spec- 
troscopic histograms now seem to be in reasonable 
agreement, though there are possible deficiencies of 
photometric objects at z ~ 0.7 and in particular, 
z-1.9. 



• Color-redshift diagrams, for the optical and NIR, 
show the close agreement between our matched 
catalo g and the "average" models of IHewett et al.l 
(|2006l ). at redshift z < 2.0. At higher redshifts, the 
models generally appear to be generally bluer than 
the mean observed quasar colors. We tentatively 
suggest this is the same affect as has rece ntly been 
reported in iWorseck fc Prochaskal (|2010f ) , namely, 
that the original SDSS color-selection that was 
used to select quasars (and ultimately our own 
matched catalog) may be systematically biased to- 
wards missing blue (u — g) < 2.0 quasars at z ~ 
2.2 — 3.5, due to the selection preferentially select- 
ing intervening H I Lyman limit systems. 

• Since stellar contamination is of great interest for 
ongoing quasar surveys such as the SDSS-IILBOSS 
and AUS, we plot a test set of stellar data with our 
matched catalog data in g JK and giK color space. 
We co nfirm findings from pr evious studies, in par- 
ticular iMaddOTeLaD (|2008l) . that (a) the stellar 
locus traces out a clear band in gJK color-space, 
related to the different stellar spectral type and 
(b that in general, quasars, especially those with 
2.2 < z < 3.2 mostly lie in a distinct region of 
gJK color space than stars. 

• We plot the iYJ colors of our matched cat- 
alog for high, z > 4.6, redshift spectroscopi- 
cally confirmed quasars, again comparing to model 
tracks, though no obvious trends are seen. Fi- 
nally, just using the ULAS DR3 and the very high, 
z > 5. 7, redshift object s reported in iFan et al l 
(120061). IJiang et all (l2008l) . IJiang et al.l (120091) and 
IMortlock et al.l (|2009[) . we find 6 (5) out of 13 (12) 
quasars have NIR (K) band detections. 



It is worthwhile to mention that ma ny of our te s ts hav e 
also been performed very recently in IWu fc Jial (|2010|) . 
and although we have not performed any direct compar- 
isons between the results found herein and their study, 
we see very good general agreement, e.g. for color-z re- 
lations, between the two investigations. 

The major motivation for the construction of the 
matched catalog presented here, was its utilization in a 
future study where we measure the if-band quasar lu- 
minosity function. As quasars are measured further in 
the NIR, the flux due to host galaxies is no longer neg- 
ligible but will rather constitute a sizeable percentage of 
the total bolometric flux from a quasar. Our future work 
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(Peth, Ross et al. in prep.), will address this, and other 
issues in order to construct and measure the observed 
-fT-band quasar luminosity function. 

Advancement can come from technological break- 
throughs as well as new theoretical insights. Future tele- 
scopes and surveys, e.g. VISTA- VHS and the SASIR, 
will not only cover more of the sky but should also be 
able to observe to greater depths. Future observations 
will warrant an updated analysis of quasar properties in 
multiple bands, and will be both necessary and essen- 
tial to further understand the formation and evolution 
of quasars. 
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APPENDIX 

A. PHOTOMETRIC BANDS AND CONVERSIONS 

Due to the differing normalizations between the SDSS and UKIDSS photometric systems, certain corrections are 
required. To present our data in the purest sense, all the NIR m a gnitu des from UKIDSS (originally AB magnitudes) 
were corrected to Vega magnitudes as suggested in lHewett et al.l (|2006l ). 

Although ULAS magnitudes are reported in terms of Vega and SDSS magnitudes are reported in AB terms for the 
most part whenever an optical-NIR color was calculated both magnitudes were left in their default term. 



Table 9 

Adapted from Table 19 of lStougfaton et al.l (f2002h . Table 1 of lLawrence et al.l p007f ) and Table 7 of lHewett eTaTI p006f) . 



O O O O 

Band A cff A A min A A max A FWHM A AB - Vega Transformations 



u 


3551 


3005 


4000 


581 


u = 


u AB - 0.927 


9 


4686 


3720 


5680 


1262 


9 = 


g AB + 0.103 


r 


6166 


5370 


7120 


1149 


r = 


t ab - 0.146 


i 


7480 


6770 


8380 


1237 


i = 


!ab - 0.366 


z 


8932 


8000 


10620 


994 


z = 


zab - 0.533 


Y 


10305 


9790 


10810 


1020 


Y = 


= Y AB - 0.634 


J 


12483 


11690 


13280 


1590 


J = 


Jab- 0.938 


H 


16313 


14920 


17840 


2920 


H = 


= H AB - 1-379 


K 


22010 


20290 


23800 


3510 


K = 


= K AB - 1-9 
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